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ABSTRACT 


Measured  elevations  in  the  wave  train  of  a  Rankine  ovoid  are 
compared  with  computed  elevations  for  an  ideal  source  and  sink. 
The  measured  elevations  are  in  phase  with  the  computed  elevations 
but  have  smaller  amplitudes.  The  discrepancy  in  amplitude  is 
attributed  to  turbulent  boundary  layer  and  wake  in  the  flow 
around  the  real  ovoid. 


ZUSAlftfENFASSPNC 

Die  gemessenen  Hohen  im  Wellenzug  eines  Rankineschen  eifcrmigen 
Korpers  werden  mit  den  berechneten  Hohen  fur  eine  i deale  Quelle  und 
Senke  verglichen.  Die  gemessenen  Hohen  sind  gleichphasig  mit  den 
berechneten,  haben  aber  kleinere  Amplituden.  Die  TJhstimmigkeit 
hinsichtlich  der  Amplitude  wird  auf  die  turbulent e  Grenzschicht  und 
die  Nachlaufstromung  am  wahren  eiformigen  Korper  zioruckge f uhrt . 

i 

RESUME 

On  compare  les  dldvations  mesurdes  du  train  dfondes  d*un  ovoide 
de  Rankine  avec  les  dldvations  calculdes  pour  une  source  et  puits 
id€ale.  Les  dldvations  mesurdes  sont  en  phase  avec  les  dilations 
calculdes,  mais  prdsentent  des  amplitudes  plus  rdduites.  La  dif¬ 
ference  d’amplitude  est  attribute  h  la  coucfae  limite  et  au  sillage 
turbulents  du  eourant  qui  entoure  1' ovoide  rdel. 
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FORTTORD 


The  material  in  this  report  was  prepared  as  part  of  the 
foundational  research  program  of  the  Naval  Weapons  Laboratory, 
Project  NO.  R360FR102,- 2101/R0110101-  Programming  of  machine 
computations  on  NCRC  was  prepared  by  Mrs.  S.  J.  Hershey. 
Analysis  and  computations  were  completed  by  7  June  1966.  The 
figures  were  prepared  or.  a  dot  plotting  cathode  ray  printer  at 
the  Naval  Weapons  Laboratory.  The  Japanese  abstract  and  the 
schematic  diagram  were  prepared  on  a  vector  plotting  cathode 
ray  printer  at  the  David  Taylor  Model  Basin. 


APPROVED  FOR  RELEASE: 


/s/  BERNARD  SMITH 
Technical  Director 
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INTRODUCTION 


I  A  project  at  the  Naval  Weapons  Laboratory  has  as  its  ultimate 

objective  the  computation  of  flow  around  surface  ships  of  finite 
?  breadth  in  water  of  finite  depth.  In  the  pursuit  of  this  project 

i  there  has  been  prepared  a  system  of  four  computing  routines,  which 

|  compute  the  components  of  velocity  in  the  wave  train  of  a  point 

source. 

A  project  at  the  David  Taylor  Model  Basin  has  had  as  its  objec¬ 
tive  the  measurement  of  wave  elevations  in  the  wave  train  of  a 
Rankine  ovoid.  When  deeply  submerged,  the  Rankine  ovoid  is  that 
streamline  of  zero  stream  function  which  would  be  generated  by  a 
simple  source  and  sink.  A  representative  streamline  is  illustrated 
in  Appendix  A. 

The  wave  elevations  as  measured  at  the  David  Taylor  Model  Basin 
present  an  ideal  opportunity  for  a  comparison  with  wave  elevations 
as  computed  at  the  Naval  Weapons  Laboratory.  One  objective  of  the 
present  report  is  to  explore  the  possibility  that  the  measured  wave 
elevations  might  indicate  how  a  sourcewise  representation  should  be 
modified  in  order  to  allow  for  boundary  layer  and  turbulent  wake. 
Before  any  conclusions  would  be  possible  it  has  been  necessary  to 
evaluate  in  detail  all  sources  of  discrepancy  between  measured  data 
and  computed  data. 

The  components  of  velocity  in  the  wave  train  of  a  point  source 
may  be  derived  from  the  Havelock  integral,  which  is  a  double  Fourier 
integral  in  wave  number  space.  The  evaluation  of  the  integral  is 
achieved  by  radial  and  azimuthal  integration  in  the  wave  number 
space.  Following  a  suggestion  from  the  David  Taylor  Model  Basin, 
DiDonato1  in  1958  developed  a  method  of  integration  in  which  the 
path  of  integration  is  so  displaced  in  the  complex  plane  as  to 
avoid  a  singularity  in  the  integrand.  Radial  integration  was  per¬ 
formed  by  Simpson’s  rule  and  azimuthal  integration  was  completed 
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with  Gauss*  rule.  Without  reference  to  this  prior  work,  Yim2  in 
1963  has  reexamined  the  problem  of  integration  in  the  near  field. 
With  the  aid  of  formulae  of  Barakat6,  or  after  radial  integration 
with  Laguerre-Gauss  quadrature,  the  azimuthal  integration  was 
completed  with  Legendre-Gauss  quadrature.  Yim  examined  also  the 
problem  of  integration  in  the  far  field  both  along  the  center  line 
behind  the  source  and  along  the  critical  line  of  wave  crests. 
However,  the  derivations  were  extended  only  as  far  as  the  first 
two  terms  of  the  asymptotic  expansion  in  each  case. 

In  the  meantime,  the  present  writer6  has  made  an  extensive  in¬ 
vestigation  of  possible  methods  of  integration.  At  an  early  date 
it  was  realized  that  radial  integration  generates  the  complex  expo¬ 
nential  integral.  In  the  sense  that  this  function  is  an  elementary 
function,  the  problem  of  integration  is  reduced  to  the  problem  of 
azimuthal  integration  alone.  The  high  accuracy  rule  for  the  inte¬ 
gration  of  a  cyclical  function  is  just  the  trapezoidal  rule  for 
equally  spaced  angles,  and  the  above  use  of  Gauss  quadrature  is 
not  appropriate-  Even  the  trapezoidal  rule  is  inefficient  for 
applications  of  importance  to  surface  ships,  and  a  new  method  of 
integration  by  parts  was  found  in  which  integration  could  be  per¬ 
formed  through  many  oscillations  of  the  integrand.  It  was  recog¬ 
nized  that  the  integrand  is  the  product  of  a  monotonic  factor  and 
an  oscillatory  factor.  When  the  monotonic  factor  was  approximated 
in  terms  of  positive  integral  powers  of  the  argument  of  the 
oscillatory  factor,  the  integration  could  be  completed  by  recourse 
to  recurrence  relations.  Further  improvement  in  efficiency  was 
achieved  through  the  expansion  of  the  monotonic  factor  in  terms 
of  positive  and  negative  half  integral  powers  of  the  argument  of 
the  oscillatory  factor.  Finally  the  asymptotic  series  for  sta¬ 
tionary  phase  with  quadratic  approximation  near  the  center  line 
or  with  cubic  approximation  near  the  critical  line  has  been 
carried  to  the  ultimate  term  of  smallest  magnitude.  These  methods 
of  integration  are  the  basis  for  the  four  computing  programs 


for  the  computation  of  wave  trains  for  a  point  source.  The  four 
programs  provide  computations  with  maximum  efficiency  or  with  con¬ 
trolled  accuracy  and  with  no  limitation  as  to  direction  or  distance 
from  the  point  source.  The  programs  were  first  reported  in  196?. 

All  of  the  above  work  is  concerned  with  the  steady  state  wave 
pattern.  In  the  integration  of  the  double  Fourier  integral  the 
classical  approach  was  to  use  the  Cauchy  principal  value  but  this 
led  to  a  train  of  waves  in  front  of  the  source  as  well  as  behind 
the  source.  Inasmuch  as  wave  trains  in  front  of  the  source  are 
unacceptable  on  the  basis  of  physical  observation,  the  singularity 
in  the  integrand  was  avoided  through  the  introduction  of  an  imag¬ 
inary  term  which  is  the  so-called  frictional  term  but  is  in 
reality  merely  a  mathematical  device  to  compel  the  path  of  inte¬ 
gration  in  the  complex  plane  to  pass  on  the  correct  side  of  the 
singularity.  An  alternative  approach  is  to  recognize  that  the 
integrand  at  the  point  of  singularity  satisfies  identically  both 
the  boundary  equation  and  the  Laplace  equation  and  can  be  added 
therefore  in  arbitrary  amount.  When  it  is  added  in  the  correct 
amount  to  eliminate  waves  in  front  of  the  source,  the  integrand 
becomes  analytic  and  the  path  of  integration  can  be  varied.  A 


more  convincing  determination  of  the  term  to  be  added  would  be 
derived  from  an  analysis  of  the  limit  approached  by  the  transient 
disturbance  which  arises  when  a  point  source  is  started  from  rest 
and  moves  with  constant  velocity  thereafter.  Starting  with  a 
formula  of  Maruo7,  an  analysis  was  made  by  Hsu8  in  1965,  who  con¬ 
cluded  that  there  was  a  transient  system  of  circular  waves  which 
were  centered  at  the  starting  point  and  a  stationary  train  of 
waves  which  trailed  halfway  back  from  the  source  to  the  starting 
point.  In  view  of  the  importance  of  the  transient  waves  as  a 
guide  to  the  properties  of  the  stationary  waves,  the  formula  of 
Hsu  is  rederived  herewith  and  is  transformed  further  to  reveal  the 
character  of  the  wave  train. 
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The  analysis  of  transient  waves  has  been  used  by  Hsu  and  Yim4  in 
1966  to  interpret  the  experimental  data  on  the  wave  profiles  of  the 
Rankine  ovoid.  Only  the  wave  profiles  on  or  near  the  center  line  were 
considered  for  interpretation  whereas  the  wave  profiles  farther  away 
from  the  center  line  are  more  interesting  although  more  unstable. 

The  experimental  wave  profiles  were  shifted  with  respect  to  time  so 
as  to  match  the  theoretical  wave  profiles  at  a  point  of  zero  elevation. 
It  was  concluded  from  a  comparison  of  the  matched  profiles  that  the 
ideal  wave  theory  predicts  wave  profiles  for  fully  submerged  bodies 
with  sufficient  accuracy  for  practical  purposes.  On  the  other  hand, 
a  systematic  discrepancy  in  amplitude  of  the  measured  ana  computed 
wave  profiles  was  noted  and  was  thought  to  be  the  result  of  viscous 
effects. 


A  series  of  determinations  of  wave  profile  have  been  made  by 
D.  A.  Shaffer8,10  at  the  David  Taylor  Model  Basin.  A.  Rankine  o.-oid 
was  towed  at  various  speeds,  and  wave  elevations  were  recorded  on  a 
Sanborn  recorder  for  various  distances  off  the  line  of  travel.  The 
position  of  the  ovoid  was  determined  from  the  angle  of  rotation  of 
a  drum  around  which  the  towing  cable  was  wrapped.  Contacts  on  the 
drum  transmitted  signals  which  were  recorded  on  the  Sanborn  record. 
Correction  was  made  for  slippage  of  the  cable  on  the  drum.  The 
model  was  provided  with  a  projecting  lug  which  tripped  a  switch  of 
known  position  and  recorded  a  step  on  the  Sanborn  record.  With  this 
calibration  the  zero  point  of  time  when  the  center  of  the  ovoid  was 
opposite  each  recording  station  could  be  computed  and  marked  on  the 
record. 


The  ovoid  had  a  length  of  4|  ft  with  a  ratio  of  length  to 
diameter  of  7.  It  was  towed  by  a  cable  of  -|n  diameter  at  a  sub¬ 
mergence  of  ly  ft  in  water  with  a  depth  of  20  ft.  The  depth  of 
submergence  was  uncertain  to  ±1. 5n  because  of  sag  in  the  cable. 
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The  length  of  travel  at  constant  velocity  was  205  ft.  In  one 
series  of  runs  the  wave  elevation  was  measured  at  0  ft  and  at 
22.75  ft  to  one  side  of  the  line  of  travel  while  in  another 
series  of  runs  the  wave  elevation  was  measured  at  11.375  ft  and 
at  34.125  ft  away  from  the  line  of  travel.  The  speeds  of  the 
ovoid  were  6,  7.3,  9,  and  10  (ft)/ (sec). 

The  amplitude  of  wave  elevation  was  observed  to  fluctuate 
with  speed  with  a  minimum  near  4  (ft)/ (sec)  and  with  a  maximum 
near  7.3  (ft)/(sec).  The  minimum  was  interpreted  in  terms  of 
interference  while  the  maximum  was  correlated  with  critical 
Froude  number. 

Seme  deviation  from  a  steady  state  wave  profile  occurred  down¬ 
stream  from  the  model  because  there  had  not  been  enough  time  for 
the  waves  to  develop  their  full  height  and  because  there  was  inter¬ 
ference  from  waves  which  had  been  reflected  off  the  walls  of  the 
basin.  The  unsteady  region  could  be  recognized  by  comparisons 
between  records  for  different  distances  of  travel.  The  length  of 
the  fully  developed  wave  train  was  reported  to  be  half  of  the 
distance  of  travel. 


STEADY  STATS 

Let  a  point  source  be  moving  at  a  depth  h  below  the  surface  with 
a  speed  0.  A  point  in  the  fluid  is  located  at  Cartesian  coordinates 
x,  y .  z  in  a  right  handed  coordinate  system  with  origin  at-  the  sur¬ 
face  over  the  source.  The  x-coordinate  is  the  distance  forward, 
tne  y-coorainate  is  the  distance  tc  the  right;  and  the  z-coordinate 
is  the  distance  downward-  The  velocity  in  the  fluid  is  expressed 


W* ( ’ 


as  the  negative  gradient  of  a  velocity  potential  <p.  The  computing 
routines  give  the  Cartesian  components  u,  v,  u  of  velocity  in  accord¬ 
ance  with  the  equations 


d<p 

u  -  -  — 
dx 


3m 

v  - - 

dy 


dq> 

v  —  -  — 
dz 


The  velocity  potential  is  the  sum  of  three  peter ti ale  in  accordance 
with  the  equation 

<P  =  +  <P*  +  %  (2) 

where  is  the  potential  of  the  source  in  an  unbounded  fluid,  c2 
is  the  potential  of  an  image  source  over  the  free  surface,  and  % 
is  a  Fourier  integral-  The  potentials  Oj  +  «p2  cancel  at  the  free 
surface,  while  the  potential  q?3  is  the  sum  of  a  monotonic  term  and 
an  oscillatory  term.  The  monotonic  term  is  responsible  for  a 
solitary  wave  or  Bernoulli  hump  over  the  source  while  the  oscillatory 
term  is  responsible  for  a  wave  train  or  Kelvin  wake  behind  the  source. 
Limiting  formulae  for  the  wave  elevation  on  the  line  directly  astern 
of  a  unit  source  may  be  derived  from  the  asymptotic  expansion  of  the 
potentials.  For  deep  submergence  the  elevation  in  the  solitary  wave 
tends  to  vary  like 


?  (x2  +  h2)s 

Fer  down  stream  the  elevation  in  the  wave  train  tends  to  vary  like 


i 

27!  2  __  k  ,  H  . 

4Xo  —  e  xor*cos(x0r - ) 

gx  4 


where  the  critical  wave  number  Xq  is  defined  by  the  equation 

_  S 

*°  0s 

The  asymptotic  formula  for  the  wave  train  is  not  useful  over  the 
source  where  it  goes  to  infinity- 
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It  is  of  interest  to  note  that  when  the  source  strength  is 
proportional  to  0  the  amplitude  in  the  wave  train  has  a  maximum  with 
respect  to  0  where 

XqA  =  |  (6) 

This  corresponds  to  a  critical  speed  of  9.8  (ft)/(sec)  when  the  depth 
of  submergence  is  1.5  ft. 

The  wave  length  X  of  the  waves  directly  behind  the  source  is  given 
by  the  equation 

.  2n 

X  =  ~  (7) 

*o 

Interference  between  a  source  and  sink  will  occur  if  their  distance 
apart  is  an  integral  multiple  of  wave  length,  while  reinforcement  will 
occur  if  their  distance  apart  is  half  an  odd  integral  multiple  of  wave 
length.  Representative  data  for  the  experimental  ovoid  are  given  in 
the  following  table 

Source  Separation  Phase  Correlation  Speed  (ft)/ (sec) 


1 

lx 

Reinforcement 

3. 7760 

E 

i 

X 

Interference 

4. 6247 

& 

I 

I 

*X 

Reinforcement 

6.5403 

It  may  be  noted  that  reinforcement  prevails  for  the  experimental  speeds. 
Computed  wave  profiles  have  been  obtained  for  the  first  two  speeds  in 
the  above  table  together  with  the  experimental  speeds.  The  speeds,  the 
critical  wave  numbers,  and  the  wave  lengths  for  which  there  are  com¬ 
puted  profiles  are  summarized  in  Table  I. 


SOURCEWISE  REPRESENTATION 


The  streamline  for  a  deeply  submerged  source  and  sink  was  deter¬ 
mined  by  iteration.  Let  the  source  and  sink  be  of  strength  ±<7  at 
the  positions  ±0  on  the  x-axis.  If  the  free  stream  velocity  is 
unity,  then  the  flux  Q  through  any  circle  of  radius  r  at  the  posi¬ 
tion  x  is  given  by  th°  equation 


_  9  1  f  x  -  a  1  if  x  +  a  } 

Q  ~  ~  ~  1  1 - zrr=  \  1 - - — - - \ 

^  l  yf  (x  -  a)2  +  r2]  ^  l  (x  4.  a)2  +  j*  J 


(3) 


The  equation  of  the  streamline  is  Q  =  0.  In  the  limiting  case 
where  r  —  0,  the  equation  of  the  streamline  is  reduced  to 


axq 

(x2  -  a2)2  " 


where  x  is  the  half  length  of  the  ovoid.  In  the  limiting  case 
where  x  —  0,  the  equation  of  the  streamline  is  reduced  to 


-a—  =  »» 

✓  a2  +  r2 


(  r  =  .1  )  (10) 

14 


where  r  is  one  seventh  of  the  half  length  of  the  ovoid.  Iterative 
solution  of  the  simultaneous  equations  (9)  and  (10)  for  a  length 
l  =4.5  ft  leads  to  the  equations 

o  =  2.0884  ft  q  -  0-3284  (ft)2  (ll) 


The  radius  r  at  any  position  x  on  the  streamline  is  determined  by  the 
Newton-Rapnson  iteration 


Q 


dr 


(12) 
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In  Table  II  are  given  recomputed  values  of  radii  for  the  same  positions 
which  were  used  in  the  specification  of  the  model. 8  The  recomputed 
values  do  not  quite  agree  with  the  original  values,  but  the  discrepan¬ 
cies  between  values  are  not  much  more  than  a  hundredth  of  an  inch  in 
the  worst  cases. 

The  partial  derivatives  of  Q  are  given  by  the  equations 


BQ 

dx 


{(x  -  a)2  +  r2}* 


{( x  4  a)2  4  r2}5 


and 


(13) 


£2  =  2ltr 
dr 


|r(x  -  a)g 
{(x  -  a)2  +  r2}1 


ir(x  4  a)q 
{(x  4  a)2  4  r2^ 


(14) 


The  components  of  unit  vectors  tangent  or  normal  to  the  streamline 
can  be  expressed  in  terms  of  the  partial  derivatives  of  Q.  Along  a 
vertical  plane  of  symmetry,  the  velocity  normal  to  the  streamline  at 
a  point  is  given  by  the  expression 


BQ  BQ 


where  the  ±  sign  depends  upon  whether  the  point  is  topside  or  bottomside. 
The  sum  of  normal  components  of  velocity  for  the  source  and  sink  is  zero 
for  deep  submergence,  but  when  the  ovoid  is  near  the  surface,  there  are 
contributions  to  velocity  from  the  surface  waves.  The  resultant  velocity 
for  the  source  and  sink  does  not  quite  give  a  vanishing  normal  component 
at  the  surface  of  the  ovoid.  The  normal  component  at  the  vertical  plane 
has  been  computed  for  the  experimental  conditions.  The  worst  case  was 
for  6  (ft)/ (sec)  where  the  normal  component  was  not  more  than  1. 6%  of  the 
free  stream  velocity. 
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FREE  SURFACE 


One  of  the  errors  in  the  mathematical  mode?  is  a  linearisa¬ 
tion  of  the  boundary  condition  at  the  free  surface.  Let  the 
equation  of  the  free  surface  be 


2  -  /(*.  y)  =  0 


(16) 


For  steady  state  conditions  the  velocity  potential  at  the  free  surface 
satisfies  the  boundary  equation 


a?  df 
dy  dy 


(17) 


and  also  the  Bernoulli  equation 


1 

2 


{< 


df  . 2  ,  df  *  ,  df 

U  +  -J)  +(-?)  +(/ 

dx  dy  dz 


(18) 


which  becomes  the  equation 


after  differentiation  with  respect  to  x.  To  the  first  order  of  small 
quantities  the  equation  of  the  surface  is 


whence  the  derivatives  of  f  are  given  by  the  equations 


I 

fi 


df  _0  d*<? 
3x~~g  JP 


9/  _  u  a2? 

dy  g  dxdy 


Substitution  of  these  expressions  of  first  order  into  the  expressions 
of  second  order  leads  to  the  second  order  boundary  equation 


a2qp  do  19  [  f  dtp  .2  ,39.2  1.99-21 

ix2  ~  *° 9z  +  1  9x  +  9y  ^  +  2  9z  j  ”° 


The  second  order  term  in  this  equation  is  the  derivative  of  a  quantity 
which  is  positive  and  fluctuates  between  minima  and  maxima  such  that 
the  minima  are  about  half  the  adjacent  maxima  in  magnitude.  The  wave 
number  of  the  fluctuation  of  the  second  order  term  is  twice  the  wave 
number  of  the  fluctuation  of  first  order  terms.  The  neglect  of  the 
second  order  term  in  the  boundary  equation  leads  to  an  error  in 
elevation  whose  amplitude  is  estimated  to  be  equal  to  the  maxima  of 
the  quantity 


g  l  dx  dy  2  9z  J 


Values  of  this  quantity  have  been  computed  for  the  experimental  condi¬ 
tions  and  the  largest  maxima  for  each  speed  are  summarized  in  Table  III. 
The  largest  error  in  elevation  is  only  about  2.0^  of  the  largest 
amplitude  of  elevation. 
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CAPILLARITY 


Another  error  in  the  mathematical  model  arises  from  a  neglect 
of  surface  tension  at  the  free  surface.  A  surface  tension  y  acts 
on  a  line  element  dr  in  the  surface  in  a  direction  which  lies  in 
the  surface  and  is  orthogonal  to  dr.  The  component  of  this  force 
in  the  direction  of  increasing  z  is  given  by  the  product  of  this 
force  and  the  rate  of  change  of  /  in  the  direction  orthogonal  to 
dr.  The  component  of  force  is  given  by  the  scalar  product 

-  yV/  •  kx  dr  (24) 

where  k  is  a  unit  vector  in  the  direction  of  increasing  z.  The 
resultant  force  on  a  closed  contour  with  position  vector  r  is  given 
by  the  circuit  integral 

-  Y  /  Vf'kXdr  (25) 

where  the  path  of  integration  is  taken  in  the  right-handed  direction 
relative  to  the  vector  k.  The  surface  tension  on  the  contour  is 
balanced  by  the  action  of  a  pressure  against  that  surface  which  is 
contain  d  within  the  contour.  Application  of  triple  product  expan¬ 
sions  and  of  the  Stokes  theorem  to  the  circuit  integral  shows  that 
the  pressure  is  given  by  the  equation 

P  =  yV  •  V/  (26) 

The  Bernoulli  equation  is  modified  therefore  through  addition  of  the 
term 


Y  ,£V 

pV3x2 
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dy2 


(27) 
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and  to  first  order  the  boundary  equation  becomes 


da<p  dtp  Y  34<p 

dP  ~  *°  Tz  "  pg  dPdP 


(26) 


The  effect  of  surface  tension  is  a  change  in  amplitude  and  wave  length 
by  the  fractional  amount 


(29) 


Inasmuch  as  the  coefficient  of  x|  in  this  expression  is  on  the  order 
of 


~  =  8  X  io~6  (ft)2  (30) 

P? 


the  effect  of  surface  tension  was  negligible  for  the  experimental 
conditions.  The  effect  of  capillarity  is  not  important  as  long  as 
the  towing  speed  substantially  exceeds  the  minimum  speed12  cf 
0.76  (ft)/(sec). 


TRANSIENT  STATE 


The  ovoid  was  started  and  stopped  in  a  tank  of  finite  size. 
Only  a  limited  part  of  the  recorded  wave  elevation  approximates 
the  limiting  profile  for  steady  state.  The  extent  of  good  profile 
can  be  deduced  from  an  analysis  of  the  transient  wave  pattern. 

Let  the  Cartesian  coordinates  x,  y,  z  of  a  point  in  the  fluid 
at  time  t  be  referred  to  a  fixed  origin  with  z  measured  downward. 
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Let  the  velocity  in  the  fluid  be  the  negative  gradient  -  Vq>  of  a 
velocity  potential  <p  which  is  a  solution  of  Laplace’s  equation 
V  *  V<p  =  0.  The  dynamical  equation  is  the  Bernoulli  equation 


at 


+ 


+ - gz  -  constant 

P 


(31) 


where  p  :s  the  density,  P  is  the  pressure,  and  g  is  the  accelera¬ 
tion  of  gravity-  Let  the  equation  of  the  free  surface  be 

*  =  /(*.  y.  t)  (32) 

and  let  a  source  be  created  below  the  origin  at  a  depth  k  and  at  the 
time  zero  in  a  fluid  initially  at  rest. 

In  the  limit  of  instantaneous  formation  of  a  finite  source  all 
terms  in  the  Bernoulli  equation  including  p/p  remain  bounded  at  the 
free  surface  and  d<p/8t  therefore  must  also  remain  bounded  at  the 
free  surface-  If  each  term  in  the  Bernoulli  equation  is  integrated 
with  respect  to  time,  then  the  term  dff/dt  becomes  the  change  in 
potential  but  the  other  terms  become  infinitesimals.  In  the  limit 
of  instantaneous  formation  of  the  source  the  potential  remains 
constant,  and  the  boundary  conditions  require  that  the  formation 
of  a  source  beneath  the  surface  be  accompanied  by  the  formation  of 
an  image  over  the  surface.  For  a  unit  source  below  the  surface  the 
potential  is  given  by  the  equation 

<Pi  =  1 . -  (33) 

✓  x2  +  y2  +  (z  -  A)2 
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and  for  an  imago  source  orer  the  surface  the  potential  q>2  is  given 
by  the  equation 


= - ==£: '  ==  (34) 

V  za  +  y2  +  (z  +  h)2 


The  Fourier  transforms  of  the  potentials  f^and  <p2  are  given  by  the 
equations 


♦*  • 
/  / 


-ic  o 


-*!*  -  *1 


♦ 


ix(%co30  ♦  >sln6) 


dude 


(35) 


and 


<p2  =  - 


_1_ 

2n 


/  / 


-it  o 


-*(*♦  fc)  ♦  iX(XC03S  ♦  >slnS) 
e  aXat7 


(36) 


To  these  must  be  added  a  potential  %  which  is  initially  zero  but 
meets  the  boundary  conditions  at  the  free  surface. 

The  linearized  Bernoulli  equation  at  the  free  surface  is 


~  +  i f  =  0  (37) 


while  the  kinematic  equation  at  the  free  surface  is 


Elimination  of  /  from  Equations  (37)  and  (38)  leads  to  the  boundary 
equation 


32<p 

at 2 


(39) 


The  Fourier  amplitude  A(x,  6 ,  t)  of  the  potential  qpa  therefore 
satisfies  the  equation 


d*A 


g 


-2  +  kca--  =  o 


(40) 


of  which  the  appropriate  solution  is 


A'x,  cf.  t)  - 


TOC 


COS  v  gX t 
rat 


(41) 


If  an  instantaneous  creation  of  a  finite  source  is  followed  after  a 
time  interval  dt  by  the  instantaneous  annihilation  of  the  finite 
source,  then  the  Fourier  amplitude  per  unit  pulse  is  merely 


A(x,  6. 


sin  ✓  gxt 


(42) 


and  from  the  boundary  equation  (37)  the  surface  elevation  per  unit 
pulse  is  given  by  the  equation 


1  r*~  —  -xfc  ♦ 

-  /  =  -  J  f  cos S  gxt  e 

71  -71  o 


ix(xcos5  ♦  ysln5 )  . 

xdxdO 


(43) 


$ 

t 

i 
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The  surface  elevation  is  initially  a  hump  but  oreaks  ultimately  into 
the  well  known  system  of  concentric  waves. 

An  alternative  formulation  is  achieved  through  the  substitution 

x  cos  G  +  y  sin  6  -  rcos  <b  (44) 

and  from  the  definition18 

1  * 

J0(z)  =  -  J  cos(z  cos  <£)a#  (45) 

71  o 

The  surface  elevation  per  unit  pulse  is  given  by  the  equation 


-  /  =  2  /  casvgKt  e  "  J0(xr)  xdii  (46) 

c 

From  the  Lipschitz  integral13  the  equation  for  the  hump  at  t  -  0  is. 
recovered.  The  initial  surface  elevation  thus  is  given  by  the 
equation 


After  a  number  of  jcnplex  transformations  cf  variable  the  elevation 
for  r  =  0  is  obtained.  The  surface  elevation  at  the  origin  is  given 
by  the  equation 


(r  =  o)  (48) 
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The  integral  in  this  equation  is  the  Fresnel  integra  of  imaginary 
argument.  The  surface  elevation  at  the  origin  subsides  gradually 
“fith  time  until  it  is  negative,  then  becomes  positive,  and  finally 


returns  to  zero. 


If  the  Bessel  function  in  Equation  (46)  is  replaced  by  its 
asymptotic  value. 


2  — 
Cs  »  A 


J0(xr)  =  ( - )2  cos  (  xr - ) 


then  there  appears  in  the  integrand  the  product 


/  K  \  - 

cos  (  xr - )  cos  V £X  t 

4 


This  product  makes  the  most  contribution  to  the  integration  when  the 
arguments  of  its  factors  differ  by  a  multiple  *  of  2n  and  furthermore 
the  derivatives  of  these  arguments  are  equal.  The  product  then 
behaves  like  the  square  of  a  cosine  and  is  positive  over  the  longest 
range  of  integration.  The  derivatives  are  equal  when 


r  =  lVl 

O  f  Y 


and  the  arguments  differ  by  a  multiple  of  2n  when 


a 

-  =  270* - 

4r  4 
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The  amplitude  of  the  waves  is  limited  by  the  exponential  function 
in  the  integrand  of  Equation  (46).  Since  neither  the  trigonometric 
function  nor  the  Bessel  function  are  ever  more  than  unity,  the 
amplitude  cannot  exceed  the  initial  elevation  at  the  origin. 

If  a  unit  source  is  created  below  the  origin  at  depth  h  and  at 
time  zero  and  then  moves  along  the  x-axis  at  constant  speed  0 ,  the 
elevation  of  the  surface  may  be  obtained  through  integration  of  the 
effects  of  pulses  which  have  been  created  during  a  succession  of 
differential  time  intervals  dx  from  T  —  0  to  T  —  t.  The  resultant 
elevation  at  time  t  thus  is  given  by  the  equation 


1  -  -*k  +  tx({X  -  0T)COS0  +  jsma) 

-  I  -  -  fff  cos  S  ?x(t  -  x)«  xdyjiBd’i 

n 


(53) 


Integration  with  respect  to  x  leads  to  a  formula  equivalent  to  that 
of  Hsu3,  while  further  development  in  terms  of  half  time  leads  to 
the  equation 


,  1  xU  cos  G  )t  -xk  +  ix{(x -J£]cose  *  jsina}  +  iSgxt 

f  =  -fj - 5__ -  e 

71  /  gx  +  xU  cos  6 


yJxde 


+  1  sin^G/'^x-  xff  cos  6  )t  g~Kh  *  ixU*  --^)eos«  ♦  yslna)  -^Sgxt 


71 


xdxd£  (54) 


v  fx.  -  xU  cos  51 


In  the  limit  as  t  —  ®  there  is  no  major  contribution  to  the  integra¬ 
tion  except  in  the  neighborhood  of  a  singularity  of  the  denominator 
of  either  integrand  where 


it  =  x  - 


cos  5 


V  ~0 


(55) 
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Local  expansions  in  terms  of  u  in  the  neighborhood  of  a  singularity 
are  given  by  the  equations 


^ ex  ±  xO cos  6  -  -  •|ff|cos^j«+  ... 


(56) 


and 

*x{(x  -  \Ut)cos6+  y  sin  0}  ±  -✓yx*  = 

2 

*jtn  . .  . 

“ — 2^  Ux_ Ut)cos9  +  y  sin0j  +  i{(x -^(Jt)cos  e  +  y  sin  0}u  +  ...  (57) 

COS  is 


Required  for  the  evaluation  in  the  limit  as  £  -  ®  is  the  remarkable 
integral 


sin  au  cos  bu 

u 


du 


f  0  if  |a|  <  | 6 | 
[  |  if  |a|  >  1 6 1 


(58) 


where  a  is  positive  but  b  may  be  of  either  sign  and  a,  b  are  identified 
as  follows: 


a  =  ^tfjcos &\t 


b  =  (x  -  f  £T£ )cos  o+y  sin  6  (59) 


The  nonzero  part  of  the  wave  train  is  confined  between  boundaries  where 

M  S  IH  - 


The  leading  edge  of  the  wave  train  is  defined  by  the  equation 

(x  -  tft)cos£  +  y  sin  6  =  0  (60) 
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and  the  trailing  edge  is  defined  by  the  equation 


4 


jga 

i 

* 

1 

Whether 

1 

depends 

phase. 

(x  -  hUt)cos  6  +  y  sin 6=0 


(61) 


d 

Te 


I  x  -  Ut  ^  y  sin  _ 
\  cos  6  cos2  6  J 


(62) 


or  wnere 


i 

■i 


„  sin 26 

y  sin  £?  =  -  — - — j-  cos  c  (x  -  Ut) 

1  +  s:n 6 


(63) 


Substitutions  from  Equation  (63)  into  Equations  (60),  (6l)  show  that 
the  wave  train  is  confined  to  the  region  where  x  satisfies  the 
inequality 


hut  cos 26  <  x  <  Ut 


(64) 


Along  the  center  line  where  =  0  the  trailing  edge  of  the  wave  train 
is  at  x  =  ^Ut  while  along  the  critical  line  of  wave  crests  where 


.-i  1 


6  =  sin'1  —  the  trailing  edge  is  at  x  =  ~Ut.  The  divergent  waves 


✓3 


extend  farther  back  than  the  transverse  waves,  as  might  be  expected 
from  a  simple  Huygens  construction. 
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When  Equation  (54)  is  so  rearranged  as  to  replace  x-^Ut  by 
x  -  Ut  the  result  is  the  eouation 


j  _  i  j»j»sin(  ^  fx-f  kJJ  cos  6)t  ^-xk  ♦  ix{(x-  ut)  cosff  *  yslnd} 

2<l  '  ....  •  ..  if _ n 


>•'  py.  +  yJJ  cos  6 


xdxdd 


<  ^  *  £¥•  ~  yJf  COS  6)t  -y.h  *  ix{ { x  -  Ut)  cos8  +  ysln#} 

+  — JJ - - - e  *  xdyJB 

c“  vgy.-vJJ  cos  6 


+  *'  +  X?  cos  @)t  -x.h  ♦  ix{  (x  •  )  cos#  +  >sintf} 

“  Z  ?X+ Xl'COS  £ 


xdxdO 


*  rrsin2-|(  »•'  fx  —  xZ7  cos  £•)  t  -xft  +  ix{ { x  -  ) cosa  •*■  >= i ns} 

JJ  ~  *  xdxdi 

71  /  rr  _ _ 


V  gx-  yJJ  cos  £? 


xjdxriB  (65) 


In  the  first  two  integrals  of  this  equation  there  is  no  major 
contribution  to  the  integration  except  in  the  neighborhood  of  the 
singularity  as  before  where  a  local  expansion  can  be  made  in  terms 
of  u.  The  integration  requires  the  value  of  the  integral 


.  sin  an  t. 

S  -  dll  -  -  (  a  >  o  )  (66) 

o  «  2 


where  a  is  defined  now  by  the  equation 


o  =  ItfjcosSjt 


(67) 


In  the  second  two  integrals  the  integrands  are  multiplied  by  the 
square  of  a  sine  which  is  always  positive  and  in  the  limit  as 
t  -  ®  has  the  effect  of  weighting  the  Cauchy  principal  value  by 

•  f- 
\  ; 

y 

* 

i 
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aasmsiWJU »■«***  ■  .  .  ’  ••wwwnstoai,  -w.^l ...  .. 


The  first  two  integrals  lead  to  the  singular  tern  which  must  be 
added  to  the  steady  state  solution  to  eliminate  wares  ahead  of  the 
source  iriiile  the  second  two  integrals  lead  to  the  Cauchy  principal 
waive  which  gives  waves  ahead  as  well  as  behind  the  source. 


SURFACE  ELEVATION 

The  components  of  velocity  for  the  steady  state  have  been  com¬ 
puted  on  NQRC  at  0. 5221  ft  intervals  in  x  from  52.21  ft  ahead  of 
the  line  of  wave  crests  to  104.42  ft  behind  the  line  of  wave  crests. 

•  The  data  for  w  at  12.0684  ft  from  the  midpoint  between  sources  were 
multiplied  by  the  quantity 


■7 

4«0 


(68) 


and  were  subtracted  to  obtain  the  computed  wave  elevations  for  the 
simple  source  and  sink.  The  data  for  x  were  converted  to  time 
through  division  by  U.  The  computed  ware  elevations  are  plotted 
against  tine  as  solid  lines  in  Figures  1  to  28. 

In  order  to  obtain  digital  values  of  the  experimental  results 
for  comparison  with  the  computed  results,  readings  have  been  taken 
at  0.1  sec  intervals  from  the  Sanborn  records  in  reference  10.  The 
records  for  133  ft  of  travel  have  been  selected  for  utilization  in 
the  present  analysis  because  the  data  for  this  distance  of  travel 
were  understood  to  be  the  most  reliable.  Some  discrepancies  in 
zero  point  settings  were  noted  and  were  called  to  the  attention 
of  ?*r.  Shaffer,  who  has  reanalysed  the  data  and  has  obtained  more 
accurate  results.  The  improved  data  are  plotted  as  dotted  lines 
in  Figures  9  to  24. 
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The  timing  of  the  experimental  data  is  believed  to  have  a  pre¬ 
cision  of  ±0.02  sec  but  some  remaining  inconsistency  suggests  that 
the  accuracy  may  be  less.  Thus  the  agreement  between  theory  and 
experiment  is  better  for  those  runs  where  the  distances  off  center 
line  were  0  and  22f  ft  than  it  is  where  the  distances  off  center- 
line  were  H|  ft  and  34^  ft.  The  discrepancy  is  especially  evident 
in  Figures  10,  14,  and  22  where  the  experimental  wave  profiles  lead 
the  theoretical  wave  profiles. 

DISCUSSION 

The  theoretical  and  experimental  wave  profiles  compare  quite 
favorably  with  respect  to  phase  in  the  range  of  fully  developed 
wave  train  and  there  is  correspondence  between  double  wave  crests 
at  higher  speed.  There  is  some  disagreement  as  to  amplitude, 
however,  part  of  which  may  arise  from  uncertainty  as  to  depth  of 
submergence.  That  this  may  have  a  substantial  effect  on  amplitude 
is  indicated  by  Table  IV  where  the  estimated  decrease  in  amplitude 
for  1.5"  increase  in  depth  is  listed.  The  effect  decreases  with 
increase  in  speed  and  is  relatively  unimportant  at  10  (ft)/ (sec). 

A  checkrun  for  this  speed  with  the  depth  increased  from  1.5*  to 
1.625*  confirmed  the  predicted  decrease  in  amplitude  and  indicated 
that  wave  shape  was  not  appreciably  affected. 

The  tests  were  run  in  water  of  finite  depth  and  at  the  highest 
speed  the  wave  length  of  the  transverse  waves  was  nearly  equal  to 
the  depth  of  water  in  the  tank.  A  checkrun  with  a  computing 
routine  for  finite  depth  showed  that  the  effect  was  negligible 
even  at  10  (ft)/ (sec)  speed. 

The  duration  of  good  wave  train  is  terminated  when  the  trailing 
edge  of  the  wave  train  passes  the  recording  station.  The  elapsed 
time  from  the  passage  of  the  model  to  the  passage  of  the  trailing 
edge  is  the  time  for  the  trailing  edge  to  advance  the  remaining 


half  distance  at  half  of  the  speed  of  the  model.  Times  of  termina¬ 
tion  in  seconds  are  given  in  Table  V.  These  values  are  for  the 
limiting  case  of  infinite  time  and  the  actual  durations  of  steady 
wave  train  tend  to  be  less. 

After  the  model  has  stopped  the  wave  train  tends  to  persist 
undiminished  except  at  the  tip  which  becomes  rounded  off.  The 
rate  of  propagation  of  the  tip  of  the  wave  train  is  reduced  to  group 
velocity  or  to  half  the  model  speed.  The  residual  wave  train  is 
reflected  from  the  far  wall  of  the  tank  but  should  not  perturb  the 
recorded  elevation  until  after  the  reflection  has  returned  to  the 
recording  station.  The  reflecting  wall  of  the  tank  was  more  than 
271  ft  from  the  starting  point  of  the  model,  and  the  residual  wave 
would  not  be  reflected  back  to  the  recording  station  in  less  time 
than  it  would  take  the  trailing  edge  of  the  wave  train  to  reach  the 
recording  station.  It  must  be  concluded  that  the  recorded  data  for 
133  ft  of  travel  offer  a  shorter  duration  of  good  wave  train  than 
the  data  for  187  ft  of  travel.  Whereas  the  data  for  133  ft  of 

travel  were  used  in  the  preparation  of  Figures  1  to  24,  the  data 
for  187  ft  of  travel  have  been  used  in  the  supplementary  Figures 
25  to  28.  The  better  agreement  between  measured  and  computed 
elevations  for  the  greater  distance  is  noteworthy. 

Error  in  the  boundary  conditions  of  the  mathematical  model 
cannot  contribute  more  than  a  few  per  cent  error  to  the  amplitude 
of  elevation.  The  assumption  of  ideal  flow  should  be  valid  at 
the  free  surface  which  is  several  radii  above  the  test  model,  but 
the  assumption  of  a  simple  source  and  sink  is  undoubtedly  in  error 
because  of  boundary  layer  and  wake.  During  the  towing  tests  the 
Reynolds  number  of  the  model  was  in  the  range  2  x  10B  to  4  x  10®. 

It  seems  likely  that  the  flow  over  the  model  was  fully  turbulent, 
inasmuch  as  the  model  was  preceded  by  a  long  towing  cable  which 
could  induce  turbulence  at  the  nose  of  the  model.  The  simplest 
effect  of  boundary  layer  separation  with  trailing  flow  in  the  wake 
would  be  a  reduction  in  strength  and  a  backward  displacement  of  the 
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sink.  The  redaction  in  strength  could  be  inferred  from  ihe  drag  on 
the  ovoid  but  there  do  not  seem  to  be  data  on  drag  for  Rankine 
oroide,  and  data  on  drag  for  torpedoes  could  not  be  extrapolated 
because  of  disparities  in  prismatic  coefficient.  Hie  drag  for  a 
Rankine  oroid  is  probably  at  least  as  great  as  the  drag  for  a  sphere 
of  the  same  diameter. 

The  presence  of  a  boundary  layer  has  the  effect  of  slowing  down 
the  real  flow  with  respect  to  the  ideal  flow  with  a  spreading  of  the 
streamlines  in  the  neighborhood.  The  spreading  becomes  pronounced 
toward  the  rear  where  there  is  separation  of  the  boundary  layer. 

The  decrease  of  longitudinal  flow  in  the  turbulent  wake  is  compensated 
by  a  decrease  of  convergent  flow  to  the  sink.  Hie  decrease  of  flow 
may  be  inferred  from  the  drag  or  rate  of  transfer  of  momentum  to  the 
fluid.  The  drag  D  is  given  by  the  equation 

co 

D  =  fPSCp  -  -2  Tip  0  f  urdr  (69) 

o 

where  S  is  the  projected  area,  CD  is  the  drag  coefficient,  and  the 
integration  of  longitudinal  speed  «  is  with  respect  to  distance  r 
from  the  axis  of  the  wake.  For  a  sphere  of  diameter  d  the  projected 
area  is  given  by  the  equation 

S=$nd2  (70) 
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At  the  experimental  values  of  the  Reynold’s  number  the  drag 
coefficient11  is  on  the  order  of  0.2  for  a  sphere  and  the  reduction 
of  strength  of  the  sink  should  be  at  least  10$  for  the  ovoid.  That 
it  is  closer  to  25$  can  be  inferred  easily  from  reruns  of  the  com¬ 
putations  with  the  strength  of  the  sink  reduced  by  appropriate 
amounts.  It  is  desirable  that  the  matter  be  settled  more  definitely 
with  the  aid  of  photographic  data  as  determined  from  smoke  traces 
in  air  tunnels  or  from  bubble  traces  in  water  tunnels. 

Since  the  preparation  of  this  report,  a  model  of  the  Rankine 
ovoid  has  been  made  at  the  Naval  Weapons  Laboratory  and  is  being 
tested  at  the  Notre  Dame  University.  Preliminary  photographs14 
with  smoke  traces  indicate  a  detectable  spreading  of  the  stream¬ 
lines  behind  the  ovoid.  Improvements  in  experimental  design  will  be 
necessary  before  quantitative  conclusions  can  be  drawn  from  the  ex¬ 
perimental  streamlines. 


CONCLUSION 


The  experimental  and  theoretical  wave  profiles  compare  quite 
favorably  in  view  of  experimental  difficulties.  The  principal  error 
in  the  theoretical  wave  profiles  arose  from  the  use  of  a  simple  source 
and  sink  for  an  approximation  of  the  ovoid.  The  sourcewise  approxi¬ 
mation  could  be  improved  from  an  analysis  of  experimental  data  on 
streamlines. 
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TABLE  I 


Speeds  for  which 

Ware  Profiles 

were  Computed 

0(ft)/(sec) 

xo(ft  r1 

X(ft) 

3. 7760 

2.2565 

2.7845 

4. 6247 

1.5043 

4.1768 

6.0 

0. 8933 

7.0337 

7.3 

0.6035 

10. 4112 

9.0 

0. 3970 

15.8267 

10.0 


0. 3216 


19.5373 


TABLE  II 


Offsets  for  Rankine  Ovoid 
with  Length  to  Diameter  Ratio  of  7 

x  {in)  r  (in) 


0.000 

3.85? 

2. 592 

3. 857 

5. 184 

3-854 

7.776 

3.850 

10. 386 

3.842 

12. 960 

3.830 

15. 252 

3.809 

18.144 

3.  757 

20.724 

3.642 

21.600 

3.570 

23.328 

3.316 

25.200 

2.669 

25.920 

2.182 

2(  400 

1.685 

27.000 


0.000 


TABLE  III 


Effect  of  Linearization  on  Surface  Elevation 
ff(ft)/(sec)  Error  (ft) 


3. 7760 

.0002 

4.6247 

.0006 

6.0 

.0017 

7.3 

.0018 

9.0 

.0018 

10.0 


.0017 


TABLE  IV 


Percent  Decrease  in  Amplitude  for  1.51 


0(ft)/(sec) 


%  Decrease 


3. 7760  28 

4. 6247  19 

6.0  11 

7. 2  7f 

9.0  5 

10. 0  4 


TABLE  V 


Times  of  Termination  of  Steady  Ware  Train 


Speed  { 

ft}/ (sec) 

Distance  (ft) 

6-0 

7.3 

9.0 

10.0 

time  (sec) 

48 

8 

7 

5 

5 

63 

11 

9 

7 

6 

95 

16 

13 

11 

10 

133 

22 

18 

15 

13 

187 

31 

26 

21 

19 

APPENDIX  C 

Plots  of  Elevation  E  (ft)  versus  Time  T  (sec) 


I- 

t 


NOTES 

1-  Figures  1  to  24  refer  to  133  ft  of  travel. 


2.  Figures  25  to  28  refer  to  187  ft  of  travel. 


3.  Dotted  lines  represent  measured  elevations  in  the  transient 
wave  train  of  a  Rankine  ovoid. 


4.  Solid  lines  represent  computed  elevations  in  the  steady  state 
wave  train  of  a  source  and  sink. 
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